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ABSTRACT 

Fifty observations at frequencies between 1.4 GHz and 43 GHz of the 6. 6-day 06.5-7+05.5-6 binary 
Cyg OB2 #5 using the Very Large Array over 20 years are re-examined. The aim is to determine the 
location and character of the previously detected variable radio emission. The radio emission from 
the system consists of a primary component that is associated with the binary, and a non-thermal 
source (NE), 0.8" to the NE of the binary that has been ascribed to a wind-collision region (WCR) 
between the stellar winds of the binary and that of a B-type star (Star D) to the NE . Previous studies 
have not accounted for the potential contribution of NE to the total radio emission, most especially 
in observations where the primary and NE sources are not resolved as separate sources. NE shows no 
evidence of variation in 23 epochs where it is resolved separately from the primary radio component, 
demonstrating that the variable emission arises in the primary component. Since NE is non-variable, 
the radio flux from the primary can now be well determined for the first time, most especially in 
observations that do not resolve both the primary and NE components. The variable radio emission 
from the primary component has a period of 6.7 ± 0.3 years which is described by a simple model of a 
non-thermal source orbiting within the stellar wind envelope of the binary. Such a model implies the 
presence of a third, unresolved stellar companion (Star C) orbiting the 6.6-day binary with a period of 
6.7 years and independent of Star D to the NE. The variable non-thermal emission arises from either 
a WCR between Star C and the binary system, or possibly from Star C directly. The model gives 
a mass-loss rate of 3.4 x 10 -5 M Q yr _1 for Cyg OB2 #5, unusually high for an Of supergiant and 
comparable to that of WR stars, and consistent with an unusually strong He I 1.083-/mi emission line, 
also redolent of WR stars. An examination of radial velocity observations available from the literature 
suggests reflex motion of the binary due to Star C, for which a mass of 23^4 M © is deduced. The 
natures of NE and Star D are also examined. If NE is a WCR, as suggested by other authors, then 
the required mass-loss rate is an order of magnitude higher than expected for an early B-type dwarf, 
and only just consistent with a supergiant. This raises the question of NE as a WCR, but its non- 
thermal luminosity is consistent with a WCR and a comparison of reddening between Cyg OB2 =#=5 
and Star D do not rule out an association, implying Cyg OB2 #5 is a quadruple system. Pursuing 
alternative models for NE, such as an unassociated background source, would require very challenging 
observations. 

Subject headings: stars:binaries — stars: early-type — stars:variables: other — radio continuum: stars 
— submillimetre 



1. INTRODUCTION 



Cyg OB2 #5 (V729 Cyg, BD +40 4220) is an eclipsing 
binary system consisting of two O-type supergiants or- 
biting in a 6.6-day perio d (jHalllll974l: [Leung & Schneider! 
ll978tlRauw et a.l.lll999h . This system is one of several lu- 



minous O-star systems in the Cyg OB2 asso ciation that 
shows evidence of variable ra dio emission (jPersi et alJ 
I1983L 119901 : iBieging et aDll989f l in observations gathered 
over ~ 20 yrs. The radio emission appears to have two 
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states: a low-flux state of ~ 2 mJy at 4.8 GHz where the 
spectral index is consistent with thermal emission from 
a stellar wind, and a high-flux state of ~ 8 mJy at 4.8 
GHz, where the spectral index is flatter than in the low 
state. The variation s appear to have a ~ 7-year period 
(|Miralles et al.|[T994f ) and have been attributed to vari- 
able non-thermal emiss ion from an expand i ng plasmon 
arising in the binary JBieging et all 119891 IPersi et all 
U99a iMiralles et alJH994h. 

Observations by lAbbott et alJ (|1981f ) with the VLA 
revealed two radio components: a primary component 
associated with the binary and a secondary radio source 
(hereaft er NE) ~ 0.8" to th e NE of the primary radio 
source. IMiralles et alJ (|1994T ) confirmed the existence of 
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NE and a 3 — 6 cm spect ral index of —0.5 ± 0.3 hinted 
at non-thermal emission. iContreras et al.l (|1997| ) argue 
that this emission is the result of a wind-collision region 
(WCR) between the stellar wind of the binary and that 
of a BO-2 V star (hereafter Star D ), 0.9" from the binary 
to the NE (jContreras et al.lll997l ). 

All previous analyses of the radio emission from Cyg 
OB2 #5 are based on observations from the Very Large 
Array (VLA), obtained in all the different configurations 
of the array, and hence data covering different spatial 
frequency (i.e. baseline length) ranges. Dependent on 
the spatial frequency coverage, the observations may or 
may not resolve the emission from both the primary and 
NE components. None of the previous analyses in the 
literature have attempted to take this into account. Fur- 
thermore, no attempt has been made to determine cither 
the continuum spectrum over a broad wavelength range 
or the location of the variable emission. As a result, 
not only is the nature of the radio emission from NE un- 
known, the evolution of the radio emission of the primary 
remains uncertain, which may account for the poor fits 
of m odels of the non-ther mal emission to the observables 
(e.g. iMiralles et aflll994l ). 

In this paper all VLA archive radio observations of Cyg 
OB2 #5 are re-examined to produce a consistently cal- 
ibrated data set. The analysis accounts for the changes 
in the resolution of the array and presents a consistent 
treatment of the emission from both the primary and NE 
sources in each observation. For the first time, the na- 
ture of the emission of both the primary and NE can be 
determined throughout the ~ 20 years of the observa- 
tions. 

2. OBSERVATIONS 
2.1. Very Large Array 

A total of 50 VLA observations of Cyg OB2 #5 ob- 
tained between 1983 and 2003 were extracted from the 
NRAO archive. These data comprise 30 observations at 
4.8 GHz (C band), 17 observations at 8.4 GHz (X band), 
2 observations at 14.9 GHz (U band), and a single ob- 
servation at 43.3 GHz (Q band). Examination of the 
few VLA observations at 1.4 GHz (L band) obtained in 
the low-resolution, C and D configurations were contam- 
inated by extended low-surface brightness emission and 
were not incorporated in this study. 

These data were edited and calibrated by a stan- 
dard approach using the NRAO's AIPS software package. 
Each observation was first examined to remove bad data. 
Phase-only antenna gain solutions were then established 
for all calibrators by self-calibration. Antennas show- 
ing random fluctuations in phase were examined closely, 
with further data editing as necessary. Both amplitude 
and phase calibration solutions were then derived for the 
calibrators. Again, further editing was applied as nec- 
essary. The absolute flux scale was established through 
bootstrapping the amplitude of the primary calibrator 
(either 3C286 or 3C84) to the phase calibrator, typically 
B2005+403 (=J2004+404), giving a flux uncertainty typ- 
ically determined to be ~ 3 — 5% 2 . Finally, phase and 
amplitude solutions of the phase calibrator were interpo- 
lated across the observations of Cyg OB2 #5. A final 

2 http : //www . via . nrao . edu/astro/calib/manual/ 



TABLE 1 
VLA Observations. 



Date 


MJD 


Cfg. 


Phase 


Primary NE 








Calibrator 


Component Component 








Flux 


Flux Flux 


(Y/M/D) 






(Jy) 


(mjy) (mjy) 




1.40Hz 


87/10/16 


47084 


AB 




2.98 ± 0.13 b 


87/11/09 


47108 


AB 


5.17 1 


3.02 ±0.14 


4.8GHz 


83/08/18 


45564 


A 


4.73 


6.21 ±0.19 0.94 ±0.10 


84/05/23 


45843 


C 


4.23 


8.15 ± 0.24 


84/05/27 


45847 


C 


3.02 2 


8.05 ± 0.26 


84/09/06 


45949 


D 


4.41 


6.81 ± 0.20 


84/09/15 


45958 


D 


4.45 


7.67 ±0.23 


84/09/20 


45963 


D 


4.34 


7.28 ± 0.22 


84/09/22 


45965 


D 


4.31 


7.55 ± 0.23 


84/09/24 


45967 


D 


4.27 


7.13 ± 0.21 


84/09/28 


45971 


D 


4.20 


6.63 ± 0.20 


84/10/16 


45989 


D 


2.94 2 


8.41 ± 0.25 


85/11/21 


46309 


D 


4.04 


5.37 ±0.16 


86/01/13 


46443 


D 


4.14 


4.18 ±0.13 


87/06/04 


46952 


D 


3.56 


3.75 ±0.16 


87/06/06 


47196 


D 


3.60 


3.90 ±0.27 


87/10/16 


47084 


AB 


3.10 1 


2.78 ±0.10 0.71 ±0.10 


87/11/09 


47108 


AB 


3.10 1 


2.57 ±0.10 0.72 ±0.10 


88/02/05 


47196 


B 


3.10 1 


4.15 ±0.22 


89/09/18 


47787 


BC 


2.76 


7.67 ±0.23 


92/11/02 


48928 


A 


2.92 


2.94 ±0.10 0.93 ±0.10 


92/12/18 


48974 


A 


3.01 


2.70 ±0.10 0.74 ±0.10 


92/12/19 


48975 


A 


3.02 


2.76 ±0.10 0.73 ±0.10 


93/05/01 


49108 


B 


3.06 


2.56 ±0.10 0.85 ±0.10 


93/12/27 


49348 


D 


3.02 


3.67 ±0.24 


94/04/09 


49451 


A 


2.90 


2.53 ±0.10 0.94 ±0.10 


94/04/16 


49458 


A 


2.90 


2.37 ±0.10 0.84 ±0.10 


95/04/27 


49834 


D 


3.00 


5.81 ±0.24 


95/06/10 


49878 


AD 


3.15 


5.13 ±0.15 0.89 ±0.10 


95/06/12 


49880 


AD 


3.14 


5.16 ±0.15 0.64 ±0.10 


97/01/04 


50452 


A 


3.08 


7.21 ±0.22 0.96 ±0.10 


00/06/30 


51725 


CD 


1.94 3 


3.30 ±0.10 


8.4GHz 


90/06/01 


48043 


AB 


3.05 


8.07 ±0.24 0.57 ±0.10 


91/10/03 


48532 


AB 


3.18 


9.03 ±0.27 0.40 ±0.10 


92/07/16 


48819 


D 


3.18 


6.99 ±0.21 


92/11/02 


48928 


A 


3.15 


4.27 ±0.13 0.66 ±0.10 


92/12/19 


48975 


A 


3.25 


3.99 ±0.12 0.57 ±0.10 


93/05/01 


49108 


B 


3.23 


4.11 ±0.12 0.39 ±0.10 


93/10/31 


49291 


D 


3.31 


4.30 ±0.25 


93/12/27 


49348 


D 


3.10 


4.90 ±0.30 


94/04/09 


49451 


A 


2.86 


3.48 ±0.10 0.52 ±0.10 


94/04/16 


49458 


A 


2.90 


3.39 ±0.10 0.65 ±0.10 


95/04/27 


49834 


D 


2.80 


7.60 ±0.23 


96/02/02 


50115 


BC 


2.88 


7.84 ±0.24 


96/12/28 


50445 


A 


2.91 


8.08 ±0.24 0.62 ±0.10 


98/07/17 


51011 


B 


2.61 


9.02 ±0.27 0.42 ±0.10 


00/06/30 


51725 


CD 


2.66 3 


3.82 ±0.16 


03/06/06 


52796 


A 


3.01 


7.17 ±0.22 0.36 ±0.10 


03/09/09 


52891 


A 


3.19 


6.09 ±0.18 0.37 ±0.10 


14.9GHz 


95/04/27 


49834 


D 


9.49 


9.39 ± 0.47 


00/06/30 


51725 


CD 


2.87 3 


5.27 ±0.26 


43.3GHz 


96/12/29 


50446 


A 


1.13 


12.17±0.61 



a Name of the phase-reference calibrator used if not B2005+403 
(=J2004+404). l=B2050+364; 2=B2200+420; 3=B2013-r-370. 

Single values given for the flux are the total flux of emission from both 
the primary and secondary components. 
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TABLE 2 
MERLIN Observations 



+41°18'20.0" — 



Date 


MJD 


Phase a 
Calibrator 
Flux 


Primary 
Component 
Flux 


Secondary 
Component 
Flux 


(Y/M/D) 




(Jy) 


(mjy) 


(mjy) 




1.40Hz 


96/02/05 
96/03/14 


50118 
50156 


3.18 
2.56 


4.57 ±0.14 
3.00 ±0.16 


1.43 ±0.14 
1.33 ±0.16 


4.8GHz 


96/11/14 


50401 


3.02 


5.82 ± 0.13 


_b 



a Phase-reference calibrator was B2005+403 (=J2004+404) for all 
observations. 

NE was "resolved out" . 

examination of the calibrated visibilities was made prior 
to image construction. 

The visibility data for Cyg OB2 #5 were deconvolved 
via model fitting using the SMERF patch (|Reidl I2006D 
to the difmap package (iShepherd et al.1 119951 ). This is 
a technique used widely in VLBI image reconstruction, 
rather than the well-known clean technique. An initial 
"best" model was established from imaging the 8.4-GHz 
visibilities obtained with the highest resolution, A con- 
figuration of the VLA. This was then used as the initial 
model for the remaining frequencies. 

As noted previously, the radio emission from Cyg OB2 
7^5 consists of two components: the primary source as- 
sociated with the O-star binary and NE, which is less 
bright than the primary. To establish models for both 
these sources, the model of the primary was first estab- 
lished. Phase-only self-calibration was then applied to 
further improve the antenna gain phase solutions derived 
initially from phase referencing alone. The model of NE 
was then established. The final model was determined by 
first fixing the location of NE and allowing the remain- 
ing model parameters to converge, before a final model 
fit with all parameters free to converge. Examples of the 
resulting deconvolved images at 8.4 GHz are shown in 
Fig.ffl 

At 8.4 GHz, the two components were readily resolved 
in all observations obtained with A and B configuration 
of the VLA, whereas at 4.8 GHz the two components were 
only resolved in A-configuration observations. In all of 
these observations, NE was always detected. In C and 
D configurations, the two components are not resolved 
separately at any of the observing frequencies, with only 
a single unresolved source being observed. 

The source fluxes were taken directly from the model- 
fitting parameters. The flux uncertainties at both 4.8 
GHz and 8.4 GHz were taken to be the maximum of 
either the rms image uncertainty or 3% of the source 
flux. At both 14.9 GHz and 43.3 GHz, the uncertainty 
was taken to be the maximum of either the image rms 
or 5% of the flux. A summary of the observations and 
observed fluxes is presented in Table [TJ 

2.2. MERLIN 

Three epochs of MERLIN observations were extracted 
from the MERLIN archive, with two observations (1996 
February 5 and 1996 March 14) at 1.4 GHz (L band) 
and one observation (1996 November 14) at 4.8 GHz (C 
band) . Editing and calibration of the visibilities was per- 
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20 n 32 m 22. 500 s 22.450 s 22.400 s 22.350 s 
Right Ascension (J2000) 

Fig. 1. — Two examples of the deconvolved VLA images at 
8.4 GHz that show the primary and NE sources. The top 
image is from 1992 December 19 during a low emission state 
and the lower image is from 1996 December 28 during a high 
emission state. From these two images it is seen that the flux 
of the primary source increases between the low and high flux 
states, whereas the flux of NE remains constant. The position 
of the primary component is consistent with the position of 
the binary system. The first contour level is 3a, where a is the 
image rms, which is 0.044 mjy and 0.028 mjy in the top and 
bottom images respectively. The contour spacings represent 
increases by a factor of 1.5 from the 3<r level. The FWHM of 
the synthesized beam is shown in the lower left corner of the 
images. 

formed using AIPS in a similar fashion to the VLA ob- 
servations. The different performance of the individual 
telescopes of the MERLIN array was accounted for by 
applying different weights to each antenna according to 
their sensitivity 3 . As in the VLA observations, the ob- 
servations were phase referenced using B2005+403 and 
the absolute flux scale was established from 3C286. The 
visibility data for Cyg OB2 #5 were deconvolved sim- 
ilarly to the VLA observations through modelling with 
SMERF. However, no phase-only self calibration was per- 
formed during the modelling process due to the small 
number of elements in the array. 

At 1.4 and 4.8 GHz, both NE and primary emis- 
sion components are resolved individually with MER- 
LIN. However, NE was heavily resolved due to the lack 
of low spatial frequency (i.e. short spacing) coverage in 

3 http : //www . merlin . ac . uk/user _guide/0nlineMUG/ 
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Cyg 



20 n 32 m 22.550 s 22.500 s 22.450 s 22.400 s 22.350 s 

Right Ascension (J20O0) 

Fig. 2.— MERLIN image of Cyg OB2 #5 at 1.4 GHz from 
1996 March 14, showing that both components are resolved. 
A taper of 250 kA was required to recover the secondary com- 
ponent. The resulting FWHM of the synthesized beam is 
shown in the lower left corner. The first contour level is 3a, 
where a is the image rms, which is 0.20 mjy. The contour 
spacings represent increases by a factor of 1.5 from the 3a 
level. 

the array. To address this issue, a Gaussian taper with 
HWHM of 250 kA was applied to the 1.4-GHz visibilities 
to increase the weight of the low spatial frequency data. 
At 4.8 GHz, NE was "resolved out" and only the primary 
component was detected. The primary emission model 
at 4.8-GHz had a major axis of 77 mas and an axial ratio 
of 0.90. A summary of the MERLIN observations is pre- 
sented in Table [2] The flux uncertainties quoted are the 
rms image uncertainties. The deconvolved image from 
1996 March 14 is shown in Fig. rj . 

2.3. Submillimetre Observations 

Two epochs of 350-GHz (850 /im) photometric obser- 
vations of Cyg OB2 #5 were obtained during CANSERV 
time with the SCUBA bolometric receiver system on 
James Clerk Maxwell Telescope (JCMT). The first ob- 
servation was taken on 1998 July 7 (MJD 51001) with 
the flux calibration established from the peak flux of the 
planetary map of Uranus. The second observation was 
taken on 1999 June 4 (MJD 51333) with Mars as the 
flux calibrator. Sky dips were used at both epochs to de- 
termine the sky opacity. Data reduction and calibration 
were performed following a standard procedure using the 
SCUBA User Reduction Facility 4 . The observed fluxes 
at 850 /im were found to be 27±7 mJy on 1998 July 7 and 
47 ± 8 mJy on 1999 June 4. Given the sizes of the uncer- 
tainties in these values, their difference does not provide 
compelling evidence for variability. In the absence of 
additional submillimetre data to suggest otherwise, it is 
assumed that these radiometry data are consistent with 
each other. 

2.4. Imaging and the NE Companion 

An infrared if-band image of Cyg OB2 #5 was ob- 
served with the United Kingdom Infrared Telescope 
(UKIRT) on 1996 May 28 using the IRCAM3 camera 

4 http: //docs . jach.hawaii . edu/star/sun216 .htx/sun216 . html 
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Fig. 3. — Maximum-entropy reconstruction of the 2.2/im im- 
age of V729Cyg oriented to show the companion to the NE 
(Star D) and drawn on a logarithmic intensity scale to com- 
press the dynamic range. 

and 5x magnifier, giving an image scale of 0.057 arc- 
sec pixel -1 . It was operated in shift-and-add mode, with 
each observation broken up into a large number of very 
short integrations. Provided that the field contains one, 
dominant point source (the central binary in this case), 
the centroid of this source on the array is located for 
each integration and the images are shifted in real time 
to bring them into alignment. Such a system is well 
suited to searching for faint companions to bright point 
sources and two such observations, each comprising 1000 
integrations of 36-msec duration, were taken of Cyg OB2 
#5. The same setup was used to observe HD 203856 at 
comparable airmass 30 minutes later to provide a point- 
spread function (PSF) calibrator. The reconstructed im- 
age is shown in Fig. [3] 

Fig. E] shows Star D to the NE, 0.90 ± 0.03" from 
the central binary at a P. A. of 55±1 °. This is con- 
sistent with the separations from the optical images 
(0.98 ± 0.06" at P.A. 61 ± 7° from CCD images and 
0.948 ± 0.043" a t P.A^ 54 ± 4° from Hipparcos) given by 
iContreras et al.l (|1997[ ). The infrared magnitude differ- 
ence between the binary and companion was measured to 
be AK = 3.1 ±0.1. From photometric observations made 
with the Carlos Sanchez Telescope (Tenerife) in 1997, it 
is deduced that the system magnitude would have been 
K = 4.52 at the time of the imagi ng observation (pha se 
0.08 on the photometric elements of lLinder et al.l <|2009f )). 
giving K = 7.6 ± 0.1 for Star D. 

If Star D is associated with the central binary, it is 
expected to have the same reddening (Ak = 0- 7, de- 
duced from Ay = 6.4 iTorres-Dodgen et al.l 119911). and 
distance modulus (10.7. ITorres-Dodgen et al.l 119911 ) as 
Cyg OB2 #5, giving an absolute magnitude Mr- = —3.8. 
This is ~ 0.5 ma g nitude too luminous for a B0 V star 
fcf. iVacca et al.l (19961 ) with intrinsic (V — K) from 
iDucati et al.l ( 2001 )). so it is suggested here that Star 
D is more evolved than a main sequence star, but proba- 
bly not a giant. The luminosity class of Star D has been 
inferred previously from photometry, yet spectroscopy is 
necessary to determine both the stellar type and lumi- 
nosity class. 

3. ANALYSIS 

3.1. Variations in Radio Emission 

The fluxes of Cyg OB2 #5 at both 4.8 GHz and 
8.4 GHz as a function of time are shown in Fig. |4j 
Through the 20 years of observations it is evident that 
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Fig. 4.— The fluxes of Cyg OB2 #5 between 1983 and 2003 at 
8.4 GHz (top) and 4.8 GHz (bottom). Fluxes from the primary and 
NE are shown as solid and open symbols respectively. Those ob- 
servations where the two components were not resolved as separate 
sources (i.e. primary+NE) are shown as triangles. The mean fluxes 
of NE at 4.8 GHz (0.83±0.11 mjy) and 8.4 GHz (0.50±0.12 mjy) 
are shown (dashed line). The MERLIN observation at 4.8 GHz is 
shown as a square. 

the 4.8-GHz emission from Cyg OB2 #5 has cycled 
through three cycles of high and low emission with an ap- 
proxim ate period of ~ 7 yrs, as first noted bv lPersi et al.1 
(1985). It is also clear that the radio emission from NE 
shows no variation and the primary radio component is 
the source of the variations in Cyg OB2 #5. 

At both 4.8 GHz and 8.4 GHz there were 12 and 
11 epochs of VLA observations respectively where NE 
was resolved separately from the primary and a flux 
could be measured directly. These observations have a 
weighted-mean flux at 4.8 GHz of 0.83 ± 0.11 mjy and 
0.50 ± 0.12 mjy at 8.4 GHz. The uncertainties are ap- 
proximately equal to our minimum image uncertainty of 
0.1 mjy and all the fluxes are within 1.7a of these means. 
Hence, the 4.8 GHz and 8.4-GHz fluxes from NE are 
taken to be constant. 

In some observations the primary and NE are not re- 
solved separately and hence only a total flux for the sys- 
tem is determined. For these observations, the primary 
source flux was determined, for the first time, by sub- 
tracting the derived mean flux observed for NE from the 
total observed flux at each frequency. 

Using fluxes from the primary at all observing epochs, 
the period of vari ation was determ ined using a string- 
length technique (jDworetskyl 119831 ). This method was 
modified slightl y to account for rela tive uncertainties in 
the flux values (|van Loo et al.|[200^ ). The period and its 
uncertainty wer e estim a ted fr om the string-length data 
as described by iFerniel (|1989f) . The string lengths as a 




Fig. 5. — String lengths for 17 8.4-GHz observations (top) 
and 30 4.8-GHz observations (bottom). The derived period 
is 7.1 ± 0.5 years for the 8.4 GHz observations and 6.6 ± 0.3 
years for the 4.8 GHz observations. Combined, these suggest 
a period of 6.7 ± 0.3 years. 

function of period are presented in Fig.[5]for both the 4.8- 
GHz and 8.4-GHz observations, along with the best-fit 
parabolas from which the periods and associated uncer- 
tainties were derived. The derived periods of the varia- 
tions are 6.6 ± 0.3 years at 4.8 GHz and 7.1 ± 0.5 years 
at 8.4 GHz. These are statistically consistent with each 
other and averaging these two results leads to an esti- 
mated emission variation period of 6.7 ± 0.3 years, which 
is remarkably similar to the "eyeball" period derived in 
previous works. 

The string-length technique was also used to search for 
periodicity in the 4.8 GHz and 8.4 GHz observations con- 
sistent with the 6.6-day orbital period of the binary. In 
this case no clear minimum was discernible. This is not 
surprising as the sample rate of the observations is much 
lower than the Nyquist frequency of a 6.6-day period. 

3.2. Spectral Indices of the Emission 

The spectral index of a continuum spectrum is a char- 
acteristic of the underlying radio emission. Assuming 
the flux S v at frequency v has a power-law behavior, 
S v oc v a , a weighted least-square fit is used to deter- 
mine the spectral index a of each component in Cyg OB2 
#5. MERLIN observations resolved the two radio com- 
ponents separately, giving a 1.4-GHz flux for NE. Since 
the two 1.4-GHz fluxes are the same within uncertainties, 
and only ~ 1 month apart, a mean flux of 1.38±0.20 mjy 
was estimated for NE. In combination with the mean 
fluxes determined at 4.8 GHz and 8.4 GHz from the VLA 
observations, a spectral index of — 0.50±0.11 was derived 
for NE (Fig. [6]). Higher frequency detections of NE would 
be useful to demonstrate if the power-law spectrum ex- 
tends to higher frequencies. 

Using both the VLA and MERLIN observations it was 
possible to determine the spectral index of the primary 
source at four different epochs of the flux variation cycle. 
Using only one MERLIN observation at 1.4 GHz and one 
VLA observation at 8.4 GHz from 1996 February during 
a high emission state, a spectral index of 0.26 ± 0.04 
was calculated. Fluxes from the VLA at 4.8, 8.4 and 
43.3 GHz during the same high flux state around 1996 
December 29 give a spectral index of 0.24 ± 0.01. In 
a fashion, the spectral index during two different low 
states was determined, with a spectral index of 0.51 ± 
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Fig. 6. — The spectrum of the secondary source NE, based 
on the fluxes of 1.38 ± 0.20 mjy, 0.83 ± 0.11 mjy, and 0.50 ± 
0.12 mjy at 1.4 4.8 and 8.4 GHz respectively. The 1.4 GHz 
observation (triangle) is from MERLIN and the 4.8 and 8.4 
GHz data (squares) are from the VLA. The best-fit spectral 
index is -0.50 ± 0.12. 
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Fig. 8. — The radio fluxes from the primary component dur- 
ing a low-emission state shown with the 350-GHz JCMT data. 
The best-fit spectrum has a slope +0.63±0.04, consistent with 
a stellar wind. The JCMT data were not obtained during the 
radio minimum but fit this spectrum. 
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Fig. 7. — The change in the continuum spectra of the 
primary radio component between the low-flux state (solid 
squares) toward the high- flux state (open squares), passing 
through intermediate states (open and solid circles). The date 
of each set of observations is marked in the upper left corner. 
The best-fit power-laws are shown for each set of data, with 
values of 0.60 ±0.04 (solid line), 0.51±0.08 (dot-dashed line), 
0.26±0.04 (dashed) and0.24±0.01 (dotted). The continuum 
spectrum during the low-flux state has a spectral index con- 
sistent with thermal emission from a stellar wind, whereas 
during the high-flux state the spectrum is flatter. It is ar- 
gued in Sec. 14. 1 1 this is due to the addition of a non-thermal 
emission component to the thermal emission from the O-star 
binary. 

0.08 on 1995 April 27 and 0.60 ± 0.04 on 2000 June 30 
from the VLA. The 15-GHz fluxes for the primary were 
determined by subtracting an extrapolated 15-GHz flux 
for NE of 0.43 mjy that assumes the power-law spectrum 
for NE holds to 15 GHz. Examination of Fig. [5] suggests 
that a 15-GHz flux around 1 mjy would present a unique 
spectrum not observed previously in a WCR. Hence, it 
is estimated that the 15-GHz flux of NE, and hence the 
primary component, has a systematic error of less than 
±0.5 mjy, which still gives low-state continuum spectra 
consistent with thermal wind emission. These data and 
spectral indices are presented in Fig. [7] The two 350- 
GHz data points are consistent with the extrapolation of 
the radio data from the low-emission state and a thermal 
stellar- wind spectrum (Fig. [8|) . 



3.3. Proper Motion and Component Separation 

A large number of observations spanning 20 years 
presents an opportunity to check for evidence of proper 
motion of the radio emission relative to the phase- 
reference source B2005±403. There is no evidence of 
proper motion in these data, with a scatter in positions 
of 100 mas, consistent with the positional uncertainty in- 
troduced by phase transfer at 4.8 GHz over ~ 5° between 
the phase calibrator and Cyg OB2 #5 using the VLA. 

Absolute position uncertainties due to phase transfer 
do not impact determination of the relative position of 
the sources in the field. Using the models from seven 
epochs of 8.4-GHz observations obtained in the VLA A- 
configuration, the position of the NE relative to the pri- 
mary was determined to be ±0.63 ± 0.07" (East) and 
±0.44 ± 0.04" (North). This implies a component sepa- 
ration of 0.77 ± 0.08" at a position angle of 55 ± 2° East 
of North, consistent with that derived from just two ob- 
servation epochs by iContreras et all (|1997| ). 

4. DISCUSSION 
4.1. The Primary and Variable Emission 

The primary emission component is associated with 
the O-star binary system and is the source of all the ob- 
served variations in the radio emission. In the low state 
the primary radio emission is found to have a spectral in- 
dex of 0.60 ± 0.04 consistent with that expected for ther- 
mal emission arising in a steady-state radially symmetric 
stellar wind. The thermal emission must be reasonably 
constant in nature since the 350-GHz observations were 
not both obtained during a radio minimum yet are con- 
sistent with the stellar wind spectrum deduced from the 
radio minimum observations, with a best-fit spectral in- 
dex 0.63±0.04 across this broad frequency range (Fig. [5]). 

For a stellar wind, the mass-loss rate can be calc ulated 
from the radio flux (e.g. I Wright fc Barlow ] [1975). As- 
suming the stellar wind has a temperature of lOkK and 
a wind composition with ionic mean charge of 1, mean 
molecular mass of 1.5 and 1 electron/ion, and assuming 
all the thermal emission arises from the binary stellar 
wind (see below for discussion of another potential source 
of thermal emission), the 4.8-GHz flux at radio minimum 
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of 2.5 mJy leads to a deduced mass-loss rate of 



M =3.4 x 10" 



1500 kms" 



1.7kpc 



3/2 



Mgyr" 



where v,^ is the term inal wind velocity, F is the vol- 
ume filling factor (e.g. lAbbott et al.1ll981h . and d is the 
distance to the source. The terminal vel o city o f 1500 
kms" 1 is adopted from iConti fe Howarth! (|1999t ) based 
on the absorption component of the P-Cygni profile of 
the He I 1.083 fim line. It should be noted that the emis- 
sion component of that line is very strong, and broader 
than those of other sup e rgian ts observed by them or in 
the atlas of iGroh et all (|2007f) . more closely resembling 
line profiles in WN-type spectra. 

The mass-loss rate derived here is consistent with 
other values determined for Cyg OB2 #5: 3.3 x 
10~ 5 M Q yr -1 from UV lines (jHowarth fc Prinjal 
1989), 2.5 x 10 -5 M yr -1 from infrared observations 
ijPersi et alj | 1990h. 3.7±1.3x Kr 5 M yr" 1 based on the 
43-GHz flux dContr eras et al"1ll996l) . adjus ted for the ter- 
minal velocity of IConti &; Howarth! (|1999| ) and for slight 
differences in adopted distance. More recently the relia- 
bility of these methods has been called into question as 
they tend to overesti mate the actual mass-loss rate if the 
wind is clumpy fe.g. lMokiem et afll2007l ; iFullerton et al.l 
2006). The most re cent estimate of m ass-loss rate in 
Cyg OB2 #5 is from Hinder et alj (|2009D . who suggest a 
rate of 2.1 ± 0.6 x 10 -5 M yr -1 for the binary system 
based on the observed rate of period change. This is in- 
dependent of distance and wind clumping but there are 
other processes which may affect the period. Neverthe- 
less, all the mass-loss rate estimates are consistent with 
each other, though notably higher than values derived 
from model atmosphere fi ts for stars of similar spectral 
type (Mo kiem et al.l 120071) . yet comparable to those of 
WR stars (ICrowtherl 120071) . This is consistent with the 
unusual strength of the Hei 1.083-/im line in Cyg OB2 
#5, indicating that the system, or one of its components, 
has a fast, heavy stellar wind comparable to those ob- 
served in WR stars. 

As the flux increases from the low state toward the high 
state the continuum spectrum flatte ns out (s ec Fig. [7]) . 
This is consistent with the findings of lPersi et al. (1990), 
who attributed the flattening to non-thermal emission 
from an expanding plasmon associated with the binary 
system. The model they produced resulted in broad 
agreement with the general variations in the data, though 
the fit to their data was poor. 

Most recently, 3D-hydrodynamical models of O-star 
binary systems with periods of a few days have demon- 
strated that radio flux variations and the flatter spectral 
index can result from variable ther mal emission a rising 
in a WCR between the binary stars (lPittardll2009f) Cv 
OB2 #5 is a contact system (|Leung fc Schneider! 1 197. 
and the nature of a wind-collision region in such a sys- 
tem is unclear. However, such a region will emit thermal 
emission that is likely variable. Whether the flux and 
variation amplitudes can be attained in such a region is 
discussed further in Sec. 14.21 

An alternative model for the primary radio emission 
component is proposed here, where the lower spectral in- 



dex during high emission is the result of the addition of 
a non-thermal component to the thermal emission from 
the binary system giving a "composite" spectrum. Such 
a model has been successfully applied to describe the rel- 
ativ ely flat continuum spe ctra of some Wolf-Rayet stars 
(e.g. IChapman et al.|[l999j ) where the non-thermal emis- 
sion arises in a WCR between the wind of the WR star 
and that of a massive companion star. 

For a system consisting of a non-thermal source em- 
bedded in a stellar wind plasma, the total observed flux 
is given as a function of frequency v and at epoch t by 

S bs(v,t) = S t h(v) + Snt(v,t) mJy. 

It is assumed the constant thermal emission component, 
S t h{v), has spectral index of +0.6 and a flux at 4.8 GHz 
of 2.5 mJy, deduced from the the primary source during 
the low emission state. Hence 



The non-thermal emission component of the total flux, 
S n t(v,t), is modelled as 



S n t(v,t) 



n 



where S^.gft) is the intrinsic 4.8-GHz flux of the non- 
thermal source at epoch t, a is the spectral index of the 
non-thermal emission assumed to be constant, and t{v 1 t) 
is the line-of-sight free-free opacity through the stellar 
wind to the non-thermal source at frequency v and epoch 
t, approximated by 



r(v,t) w T 4 . 8 (t) 



4.8/ 



where T4.g(i) is the 4.8-GHz line-of-sight free- free opacity 
at epoch t. 

The line-of-sight opacity is dependent on the geometry 
of the line-of-sight to the non-thermal emission. Here, 
the case of a n on-thermal sour c e in o rbit about the binary 
is considered. iWilliams et al.l (|1990f ) derived the varying 
free-free opacity along a line-of-sight to a non-thermal 
source orbiting in the circumbinary wind of the massive 
WR+O binary WR 1 40. Following Eqns. 12 and 14 in 
IWilliams et al.l (|1990f ). the opacity is dependent on the 
orbit inclination (i), argument of periastron (uj), as well 
as the epoch-dependent true anomaly (/) and the sepa- 
ration of the orbiting source from the companion star (r) 
in units of semi-major axis distance (a) such that 



74.8 (*) 



£seci 



2Ar 3 cos 3 (w + /) 



sin(oj + /) cos(w + /) tani+ 



(1 + tan 2 i) arctanl 



(1) 



A 



where 



\tan(w + /) tani 



1 + tan 2 (u + /) + tan 2 (i) 



and £ is a constant proportional to the square of the ion 
density in the stellar wind at a ra dius equal to the semi - 
major axis a, and related to in IWilliams et al.l (|1990f ) 
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Fig. 9. — The best-fit orbiting non-thermal source model 
is shown against the observed fluxes of the primary at 8.4 
GHz (top) and 4.8 GHz (bottom) for the s = 0, 0.5, 1, and 2 
models. 
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The intrinsic non-thermal flux 5*4.8 (t) is expected to de- 
pend on the local conditions e.g. electron density, which 
will vary as the source moves through the dense circumbi- 
nary wind. This may be approximated by assuming a 
simple power-law relation with separation, namely 

S 4 . 8 (*) =S 4 . 8 r- s , 

where S'^ 8 is the non-thermal flux when the separation 
is equal to a, and s is the power-law index. These dcfini- 
tions, along with the analytic solution to equation [T] (cf. 
Williams ct al. 199(3), allow S b s (v, t) to be determined as 
a function of the orbital phase of the non-thermal source 
orbiting the binary system. 

A standard Levenberg-Marquart x 2 -minimization 
technique was applied to both the 4.8-GHz and 8.4-GHz 
fluxes of the primary component to determine values for 
the seven free parameters in the model for each of the 
cases s = 0,0.5,1, and 2, assumi ng an orbital period of 
6.7 years. [Doughert y et alj (|2003l ) suggest s = 0.5 for the 
non-thermal luminosity of a WCR. The resulting model 
parameters are given in Table [3] and the light curves aris- 
ing from these models are plotted in Fig. [9] Fig.fTOlshows 
the models for the s = and s — 0.5 cases folded into 
the 6.7-year period. 

Each of these simple models show a good fit to the data 
and are effective in matching observations at both fre- 
quencies across all four observed emission cycles through 
the 20 years of observation. The reduced-x 2 values range 
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Fig. 10. — The best-fit orbiting non-thermal source model is 
shown for the s = (top) and s — 0.5 (bottom) cases at 4.8 
GHz (solid line) and 8.4 GHz (dashed line) against the fluxes 
of the primary, phased with a period of 6.7-year. Parameters 
for each model are given in Table [3] 



from 4.4 to 4.6 indicating that each model has its flaws, 
with the s = case being formally the best-fit. The s = 
case corresponds to intrinsic non-thermal emission that is 
constant throughout the orbit i.e. the variation is caused 
entirely by the varying line-of-sight opacity. In this case 
the best-fit model is a good fit to the observations except 
it is too low to match the peak emissions at 8.4 GHz. 
Introduction of non-zero values for s gives models able 
to match these sharp emission peaks but the resulting 
model flux is too low to match the 4.8 GHz observations 
when the primary flux is in decline at MJD 46309 and 
MJD 46443 from the high to low emission state, just prior 
to radio minimum (arou nd deduced orbital phase ~ 0.1). 

Willi ams et al.l (|1990) could not fit the radio flux varia- 
tions in WR 140 with this model and a single value of Cg 
(corresponding to our £) . This was attributed to the very 
different densities of the WR and O-star winds traversed 
by the line-of-sight at different orbital phases. The qual- 
ity of the fits to the Cyg OB2 #5 radio-flux variations 
with a single value of £ suggests that the stellar winds in 
this system have comparable densities, consistent with 
stars having comparable mass-loss rates rather than the 
~ 30-fold ratio between the WR and O-star winds in 
WR140. 

In all four models the stellar wind is completely opaque 
during the low emission state. Unfortunately the ra- 
dio observation at both frequencies are not evenly dis- 
tributed across the 20 years of observations and there 
are large gaps in the phase coverage of the observed light 
curves. Consequently the model fitting is unevenly fo- 
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TABLE 3 

Best model-fit parameters for the orbiting non-thermal source model. 



s 


, S 4.a . 
(mjy a s ) 


ct 


(°) 


i 

(°) 


e 


To 
(MJD) 








0.5 
1 
2 


5.3 ±0.5 

6.4 ±0.6 
7.1 ±0.6 
7.6 ±0.7 


— 0.18 ± 0.25 
— 0.34 ± 0.26 
-0.42 ±0.28 
-0.47 ±0.30 


319 ±3 
315 ±3 
352 ±5 
23 ±5 


90 ±40 
90 ±25 
88 ±46 
85 ±48 


0.69 ±0.04 
0.44 ±0.04 
0.23 ±0.04 
0.11 ±0.05 


53836 ± 35 
53516 ± 34 
53636 ± 30 
53737 ± 31 


0.48 ±0.13 
0.62 ±0.12 
1.36 ±0.22 
1.64 ±0.25 



cused by the group of observations obtained during the 
low emission states. Based on the s — model, the 
last low emission state occurred around 2007 February, 
and the next high emission state will occur around 2010 
November. 

The inclination across all models is consistent with 
an eclipsing orbit, though this parameter is poorly con- 
strained. Likewise, even though the non-thermal spectral 
index tends to become more negative with increasing s, it 
is also poorly constrained. The non-thermal flux S' A g in- 
creases with increasing s as would be expected: a greater 
flux fall-off with separation would require stronger emis- 
sion to match the observations. Varying s ranges the 
time of periastron passage (T ) by up to a year between 
the s = and s = 0.5 cases, as evident in Fig. QUI The 
remaining To values fall within this range. By far the 
greatest impact of varying the s parameter is on the de- 
duced orbit eccentricity. For s = 0, the orbit eccentricity 
is high with e ~ 0.7, and as s increases the eccentricity 
decreases sharply, with e ~ 0.1 for the s = 2 case. As 
noted above, many of the observations were made dur- 
ing the low-state when the stellar-wind is opaque and 
the non-thermal emission does not contribute to the ob- 
served flux being fit by the model. This contributes to 
the uncertainties in the fitted parameters. 



4.2. Evidence for a Third Star? 

A non-thermal source orbiting the binary system re- 
quires a star (hereafter Star C) to be in a 6.7-year orbit 
around the binary. This star could contribute the non- 
thermal radio emission via a WCR arising from the colli- 
sion of its own s tellar wind with the w ind from the O+O 
star binary fe.g. lEichler fc UsovHl 993). Such WCRs have 
been observed directly i n some WR+O star and O +O 
star binary systems fe.g. lDoughertv fc Pittardll2006L and 
references therein) . Alternatively, the non-thermal emis- 
sion may arise from the putative third star directly, e.g. 
a compact object. 

Given the high luminosities of the two supergiants in 
the binary and emission from circumstellar material, it 
will be very hard to detect the proposed third star di- 
rectly, let alone measure its orbit. Instead, the radial 
velocities (RVs) of the central binary are examined to 
search for reflex motion due to it being in an orbit with 
Star C, as suggested by t he radio obs e rvatio ns. 

The RVs measured by IRauw et al.l (|1999f) come from 
four observing runs, each between 1 and 4 weeks duration 
and separated by about a year. As the variations should 
coincide with a period near 6.7 years, each of these runs is 
treated as a single observation. A fift h observation comes 
from t he first five RVs measured by iBohannan fc Contil 
( 1976) in the space of a month 23 years earlier. For each 



TABLE 4 

Mean RV deviations (O-C) from the 
O+O orbit as a function of radio phase 

FOR FIVE OBSERVING RUNS. 



MJD 


4> 


n (RVs) 


mean (O-C) 


CT 


41150 


0.81 


5 


14.9 


9.0 


49567 


0.26 


11 


-7.6 


6.7 


49914 


0.40 


7 


-2.2 


6.6 


50316 


0.56 


3 


-3.7 


7.3 


50640 


0.69 


1 


8.6 


8.7 



RV observed from the primar y 5 , the resid u al (O -C) was 
calculated from the orbit bv IRauw et alJ (|1999l ) (based 
on all the RVs) and formed the average (O-C) for each 
run. These are given in Table 01 together with the radio 
orbital phases calculated using P = 6.7 years and To for 
the s = model from Table [3] 

A systematic increase of RV between phases 0.26 and 
0.89 is seen, implying that the O+O binary moves away 
from us more rapidly. This implies that Star C moves to- 
wards us more rapidly in this phase interval so that the 
circumbinary extinction to the non-thermal radio source 
diminishes, consistent with it brightening during this or- 
bital phase. 

The run of mean (O-C) with phase is compared with 
the reflex motion of the O+O binary in orbit with Star 
C following the orbital elements of the embedded non- 
thermal radio source from the s = case (see Fig. ITT]) . 
Fitting 

v r = 7 + Ko+O (e cos u + cos(/ + uj)) 

for Ko+o and systemic velocity 7, gives i^o+o = 32 ± 
17kms _1 and 7 = — 5. 9+4. 7 km s _1 , both very uncertain 
given the uncertainties in the (O-C)s and the shape of 
the RV curve in the phase range of the observations. The 
non-zero 7 is a consequence of not weighting the mean 
(O-C)s by the rather unequal numbers of observations 
from which they were deduced. The mass function f{m) 
can be derived from P (in days) and K (in kms -1 ) from 

(mo+o + mc) 

=1.036 x 10- 7 (1 -e 2 ) 3/2 iv 3 P, 

allowing an estimate of the mass, m^, of Star C. From 
the data here, f(m) = 3.2+2's M©. Assumi ng sin(i) = 1 
and adopting mo+o = 41.5 ± 3.4 M© from lLinder et all 
(|2009|) . this gives m c = 23±f 4 M Q for Star C. The large 

5 The secondary has not been used as IRauw et al.l 1119991 1 find 
a significantly lower 7-velocity for it in their orbit solution and 
attribute this to formation of the absorption lines in the wind 
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Fig. 11. — Comparison of observed (O-C) residuals, plotted 
with 2-a error bars against radio phase, and the RV curve for 
the reflex motion corresponding to the elements of the s — 
model (i.e. with u> shifted by 180° to reflect the location of 
the binary system in the orbit, rather than the 3rd star), with 
K = 32 ± 17 kms" 1 and 7 = -5.9 ± 4.7 kms" 1 giving the 
best-fit of the observables. 

uncertainty in the mass stems from the high relative un- 
certainty in K and the X 3 -dependence of f(m). 

With a paucity of observations, most especially at 
phases of the putative orbit where the radial velocity 
changes most dramatically, the uncertainties in this anal- 
ysis are high. However, until further observations can be 
obtained to test this analysis, it remains a tantalizing 
piece of evidence that Cyg OB2 #5 is a tertiary system 
rather than a binary. 

Additional support for Star C comes from X-ray data 
that reveal a hard component. iLinder et alj (|2009f ) sug- 
gest this is likely to arise in a WCR and argue that in 
such a compact binary like Cyg OB2 #5 the stellar winds 
of each component are far from achieving terminal veloc- 
ity and the resulting X-ray emission would be rather soft. 
Hence they suggest a WCR would be between the binary 
and another star. Such a WCR between the binary and 
Star C provides a ready source for the non-thermal ra- 
dio emission in the model described in Sec. 14.11 Such 
a mechanism has been widely established for WR stars 
and many O star syst ems that exhibit non-thermal emis- 
sion (|De Beckerll2007t) . A mass of ~ 23 M for Star C is 
consistent with a late O/early B-type star, which would 
have a sufficiently strong stellar wind to produce a WCR 
with the binary system wind in an orbit of size ~ 14 AU. 

The recent models of thermal emission in WCRs by 
IPittardl ([2009) raise the possibility that variable plasma 
density in a WCR between the binary and star C might 
account for the variable radio emission. Given the sepa- 
ration of the 6. 7- year orbit, a WCR between the binary 
and Star C is undoubtedly adiabatic (certainly away from 
orbit periastron) and thus any thermal free-free emission 
in the WCR will be optically thin, with a spectral index 
of —0.1. In combination with the stellar wind continuum 
from the binary system and Star C, this could result in 
a continuum spectrum that is flatter than a stellar wind 
spectrum if the thermal flux from the WCR is s ufficie ntly 
high. The fluxes in the simulations of IPittardl (|2009f ) are 
two to three orders of magnitude less than observed in 
Cyg OB2 #5, attributable to the use of mass-loss rates 
~two orders of magnitude less than deduced here for Cyg 
OB2 #5. The optically-thin thermal flux scales as the 



total number of ions in the WCR, which is oc M 2 /D, 
where D is the distance from the binary to the WCR. 
Since the deduced mass-loss rate and separation in Cyg 
OB2 #5 are respectively ~ 100 and ~ 40 times thos e 
used in the simulations (cf. model cwb2 in Pittard 2009) , 
the optically-thin flux from the WCR could be a similar 
order of magnitude as the observed fluxes. However, in 
this adiabatic scenario the flux variations scale as D^ 1 , 
and a highly eccentric orbit, as deduced in the s = 
model presented above, spends the bulk of the orbit near 
apastron where D changes little. Hence, the flux would 
change little, contrary to the observations. Thus variable 
thermal emission from the WCR alone can not account 
for the radio light curves in Cyg OB2 #5, though free- 
free opacity through the circumbinary wind undoubtedly 
plays a role. This possibility needs to be explo red fur- 
ther, though the types of models described by IPittardl 
(2009]) are beyond the scope of this paper. 

Lastly, an alternative source of the X-ray and non- 
thermal radio emission could be a compact object such 
as a neutron star, though the estimated mass from the 
reflex motion analysis implies this possibility is remote. 
Though it is not clear how a compact object produces 
the radio emission, stellar wind accretion onto a ~ 2 M Q 
compact star could generate the observed X-ray luminos- 
ity if gravitational potential can be converted to X-ray 
power efficiently. The gravitational capture radius of a 
neutron star of mass M n for a stellar wind of velocity 
Voo = 1500 kms^ 1 is given by 

= ^GMn ^ 2 3 x 1Ql0 cm 

Assuming the gravitational energy of the captured stellar 
wind is converted to X-ray emission with efficiency e, 
then luminosity would be 

GM n nB? g F m (r)e 
L x = -2 3.1 x 10 41 .F m (r)e erg s \ 

tin 

where F m (r) is the mass flux of stellar wind at distance 
r. For an orbit major axis of 14 AU and an eccentricity 
of 0.7 (s = model), periastron separation is 4.2 AU at 
which distance F m = 4.3 x 10~ 8 gcm~ 2 s _1 for a stellar 
wind mass-loss rate of 3.4 x 10~ 5 M Q yr _1 . This gives 
L x ~ 3.4 L w for e = Inconsis tent with the L x — 1.5 L Q 
derived bv lLinder et al.1 (|2009l ) and adjusted to a distance 
of 1.7 kpc. 

4.3. The Secondary Source, NE 

The spectral index of — 0.50 ± 0.11 derived here indi- 
cates the radio emission fro m NE is non-th e rmal. This 
was previously suggested bv lContreras et all (|1996[ ) from 
a limit to the spectral index of —2.4 ± 0.6 deduced from 
one e poch of observat i ons at 4.8 and 8.4 GHz in 1994 
April. iContreras et alj (|1997l ) were the first to note that 
NE lies directly between the Cyg OB2 #5 binary and a 
B-type star (Star D), ~ 0.9" to the NE. This led them 
to propose that NE is the result of a WCR between the 
stellar wind of the Cyg OB2 #5 binary system with that 
from Star D. They further argued that the separation of 
NE from both the primary and Star D was consistent 
with the expected relative wind momenta of the binary 
and the B-type star. This assertion is re-examined here 
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based on seven epochs of 8.4-GHz observations of the 
primary and NE. 

The position of the WCR relative to the positions of 
the sources of the colliding winds is given by 

(e.g. lEichler fe Usovlll993l ) where D is the separation be- 
tween the primary O-star binary (or triple) system and 
the B star, ro+o is the distance from the WCR to the 
primary system and 77 is the the wind-momentum ratio 
of the two stellar winds given by 

M D v D 

V = —■ • 

Md+o v o+o 

Here Mo+o> Mr>, vo+o an d vd are the mass- loss rates 
and terminal-wind velocities of the O-star binary and 
Star D respectively. The separation of the O-star bi- 
nary and the B-type companion is D = 0.93 ± 0.02" as 
determined from a weighted average of separations de- 
duced from optical and IR observations (Sec. 12.4( 1. Com- 
bined with ro+o = 0.77±0.08" measured here (Sec. l3.3j) . 
r\ = 0.04^q q3, where a Monte-Carlo method was used to 
determine the uncertainty in 77 as it is an ill-behaved 
function. Adopting Mq+q = 3.4 x 10~^ yr' 1 with 
v 0+ o = 1500 kms -1 (jConti fc Howarthlll999h leads to 

M D = 0.5-6.1 x 10- 6 ( 1Q0Q v ^ 1S " 1 ) M yr -1 . 

Keeping in mind the mass-loss rate of the binary is high 
compared with values deduced from model atmospheres, 
and the wind terminal speed for Star D is unknown, it 
is noted that only at the lower extremum is the mass- 
loss rate consistent with those expected for late-O/early- 
B supergiants, with early-B dwarfs ha ving mass-loss 
rates around an order of magnitude lower (Mokiem et al. 
120071) , though it is noted that the IR photometry of Star 
D suggest a star that is more luminous than anticipated 
for a B0 dwarf. The need to consider the extremum 
mass-loss rate in order to account for the relative loca- 
tion of the binary, NE and Star D raises the question of 
whether NE is truly a WCR. However, the Wolf-Rayet 
system WR 147 provides a ready example of an early B- 
type dwarf star providing a sufficiently dense wind to 
give a readily obs erved WCR with the dense wind of a 
WN8 companion (|Williams et al.lll997f ). 

The proximity (in space) of Star D and the Cyg OB2 
#5 binary can be tested by comparing their respective 
reddenings. Combining the K magnitude determined in 
Sec. 12.41 with the visual magn itude V = 13.1 ± 0.4, de- 
rived bv lContreras et all (|1997D . gives (V-K) = 5.5±0.4 
for Star D, assuming it does not vary. This implies a red- 
dening of Ay — 6.9 ± 0.5 for Star D, greater than that 
(A y = 5.7±0.3) implied b y (B-V) = 1.6 ±0.1 measured 
by IContreras et all (|1997l ). In spite of this apparent dis- 
crepancy between visual and visual/IR-determined red- 
denings, the fact that the reddening estimates bracket 
that of Cyg OB2 #5 (Ay — 6.4), there is no reason to 
rule out an association between the binary and Star D, 
and hence the possibility that Cyg OB2 #5 is a quadru- 
ple system. 



It is possible to estimate the luminosity of a WCR 
based upon the kinetic energy of the two colliding winds. 
The surface area of the WCR can b e approximated 
as a spherical cap of diameter irrn (Eichlcr fc Usovl 
119931 ) where rrj = D — ro+o is the distance from 
Star D to the WCR. Using this area combined with 
the kinetic luminosity of the stellar wind of the bi- 
nary, Lo+o = 2.4 x 10 37 ergs -1 , leads to luminosity 
Lwcr = 6.4 x 10 35 ergs" 1 entering the WCR. The ra- 
dio synchrotron luminosity L syn arising from the WCR, 
is estimated from L sim ~ 10~ s L wcr (jChen fc White! 
[1991 IPittard fc D ougherty! l200?l . giving an estimated 
synchrotron luminosity from the WCR of ~ 6 x 10 27 
ergs -1 . The radio luminosity of NE is estimated by in- 
tegrating the observed continuum spectrum. Assuming a 
power-law spectrum between 0.1 - 100 GHz gives a syn- 
chrotron luminosity of 2.5 x 10 26 ergs -1 for a distance 
of 1.7 kpc. Considering this is an order-of-magnitude ar- 
gument, the synchrotron luminosity of the NE is closely 
consistent with that anticipated from a WCR. 

An alternative model for NE is that of a background 
radio source in chance alignment with Cyg OB2 ^5 and 
Star D. It is difficult to refute this possibility unequiv- 
ocally, though the probability of such an alignment oc- 
curring randomly within 1 arcsecond of Cyg OB2 ^5 is 
very low. Extragalactic source counts at 1.4 G Hz indi- 
cate ~ 70 sources deg -2 of around 1 mJy (e.g. IJacksonl 
I2005T ). implying - 10" 6 of these sources in 1 arcsecond . 
If the source is at a much greater distance than Cyg OB2 
#5 it will have a different proper motion to the binary. 
A distant galaxy will remain fixed relative to the refer- 
ence frame of background quasars while the binary and 
Star D will exhibit proper motion relative to the frame. 
At radio wavelengths, it may be possible to measure this 
proper motion through astrometry. Certainly, the proper 
motion determined from the VLA observations discussed 
here is less than 100 mas, and to improve this precision 
would require VLBI observations. The low radio bright- 
ness of NE presents a challenge for higher precision VLBI 
astrometry. Alternatively, IR/optical imaging could re- 
veal if NE is associated directly with an object. Cyg 
OB2 #5 is bright in the optical (V=9.21) and attempts 
at imaging the region between the binary system and 
Star D are thwarted by the high contrast of the binary 
and the large PSF of the imaging telescopes. The IR im- 
age (Sec. 12. 1} was searched for evidence of a background 
source at the location of NE, but with no success. A 
smaller PSF could be attained by interferometry, but 
the contrast could only be defeated through either adap- 
tive nulling or coronographic imaging. This observing 
challenge remains to be attempted. 

5. SUMMARY 

This paper re-examines over 50 VLA observations of 
the well-known O-star binary system Cyg OB2#5 in an 
attempt to locate and characterize the well-known vari- 
able radio emission in the system. The radio emission 
consists of a primary component that is associated with 
the binary system and a component to the NE. Both 
components are resolved in all 23 epochs of highest res- 
olution VLA A-configuration observations, which reveal 
the flux of NE is constant while the flux of the primary 
varies. The constant flux from NE permits the flux of 
the primary to be derived in all observations for the first 



12 



time, most especially those where the two components 
are not resolved individually. A string-length analysis of 
the derived radio light curves of the primary emission at 
both 4.8 GHz and 8.4 GHz gives a period of 6.7±0.3 years 
for the variations. 

The primary emission changes character as it varies be- 
tween its high state of ~ 8 mJy at 4.8 GHz and a spectral 
index of 0.24±0.01 to a low emission state with a 4.8 GHz 
flux of ^2 mJy and a spectral index of 0.60 ± 0.04, con- 
sistent with thermal emission from a stellar wind. Ob- 
servations at 350 GHz obtained at an epoch not during 
a low emission state are also consistent with the thermal 
emission level observed during the low state, and hence 
lend support to the argument that the thermal emission 
observed during the low state remains relatively constant 
through the orbit. 

The mass-loss rate of the binary is deduced to be 
3.4 x 10 -5 M yr" 1 from the flux during the low emis- 
sion state, which is unusually high for an Of supergiant, 
and comparable to the rates determined for WR stars. 
Together with the anomalous strength of the Hei 1.083- 
/jm emission line, also redolent of a WR star, this points 
to a fast, heavy w i nd an d supports the suggestion by 
IBohannan fc Contil (fl976l ) that Cyg OB2 #5 is an im- 
mediate progenitor of a WR binary system in which the 
mass loss has not yet revealed enough evolved core ma- 
terial to affect the observed spectrum. 

The flatter spectral index during the high state is at- 
tributed to the addition of a non-thermal component to 
the thermal emission from the binary stellar wind. A 
non-thermal source orbiting within the stellar wind en- 
velope of the binary system every 6.7 years can account 
for the variations in radio flux through orbit modulation 
of the free-free opacity along lines of sight to the non- 
thermal source. 

Such a model requires the presence of a third star in 
association with the binary system. The high luminos- 
ity of the binary components and the emission from the 
stellar wind make the detection of a third star challeng- 
ing. An analysis of radial velocity data from the liter- 
ature provides supporting evidence of reflex motion in 
the binary as a result of a third star, labelled Star C, 
with a mass of 23l^ Mq. Until further observations 
can be obtained, especially at orbital phases where the 
radial velocities change most dramatically, this provides 
a tantalizing piece of evidence that Cyg OB2 #5 is ter- 
tiary system rather than a binary. Additional support for 
a third star comes from a hard X-ray component, that 
iLinder et all (|2009T) suggest arises in a WCR between the 
binary wind and that of a third massive star. 

This study also re-examines the NE source and its 
previous identification with a WCR between the winds 
of the O-star binary and that of a B0 star (Star D) 
0.9" to the NE. Using MERLIN observations at 1.4 GHz 
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with the VLA observations confirms the non-thermal na- 
ture of NE with a spectral index of —0.50 ± 0.12, and 
gives the relative separation of the binary and NE to be 
0.77 ± 0.08". Higher frequency radio observations would 
be useful to reveal the properties of the underlying rel- 
ativistic electron population. Through wind-momentum 
balance, the mass-loss rate of Star D is estimated to be 
between 0.5 — 6.1 x 10 -6 M Q yr _1 , consistent with a late- 
O/early B supergiant at the lower extremum of this mass- 
loss rate range, and an order of magnitude too high for 
a lower luminosity star. This raises the possibility that 
NE is not a WCR, though the WR+B binary WR 147 
provides an example of an early B-type dwarf providing 
a sufficiently dense wind to give a readily observed WCR 
with the dense wind of a WN8 companion. Analysis of 
IR observations of Star D to the NE reveal an apparent 
discrepancy between visual and visual/IR-based redden- 
ing estimates, but provide no compelling reason to rule 
out an association between Cyg OB2 #5 and Star D, 
and hence the possibility that Cyg OB2 #5 is a quadru- 
ple system. An estimate of the non-thermal luminosity 
of NE is also consistent with a WCR. To test the alter- 
native possibility that NE is an unassociated background 
source requires either high precision proper motion ob- 
servations through VLBI, or deep optical IR imaging. 
Both of these possibilities require very challenging ob- 
servations, that remain to be attempted. 
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